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ABSTRACT 


A new hybrid multilevel-linear quadratic regulator ( ML-LQR) 
approach was developed and applied to the attitude control of 
models cf the rotational dynamics of a prototype flexible 
spacecraft and of a typical space platform. Three axis rigid body- 
flexible suspension models were developed for both the spacecraft 
and the space platform utilizing augmented body methods. Models of 
the spacecraft with hybrid ML-LQR attitude control and with LQR 
attitude control were simulated and their responses with the two 
different types of control were compared. 
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SECTION 1 


L.<) INTRODUCTION 

This report s submitted in compliance with the i-cop 1 jf 
W >rX under contract NAS8-3J979, Tb« period of 
performance covered by the contract is from August 15, 
1980 to October 15, 1981, The submission and approval of 

this report constitute the successful completion of the 
"Exhibit A" portion of the contract. 

This reoort is a sequel to two others previously 
submitted under a different contract number. The two 
prior reports, reverences (1-1) and (1-2) were submitted 
in October 1978 and September, 1979 and covered the 
periods from July 27, 1977 to July 27, 1978 and r rom 
August 2b, 1978 to Adjust 2b, 1979, respectively, if', 
compliance with "Exhibit A M of contract NAS8-32bbb. 

1 JBJKCTI YES 


Too sections that follow summarize the effort expended on 
che Modular Design Attitude Control System Study 
contract. The primary objective of the study was to 
develop a new approacli to applying attitude control to 
mathemc.t ical models of the rotational dynamics of a 
representative flexible spacecraft and of a 
repress ntat ive space platform incorporating both the 
flexible spacecraft and a deployable truss. Since these 
rotational dynamics models were comprised of equations 
containing second order derivatives, it was necessary to 
develop a new control approach th*t could treat such a 
Claris of models effectively. A secondary objective was 
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the development of more efficient methods for generating 
the three axis discrete mass models of the rotational 
dynamics of the flexible spacecraft und space platform. 
These objectives led to the development and application 
of a new hybrid control approach combining the techniques 
of multilevel (ML) control and an extension of linear 
quadratic regulator (LQR) control and the application of 
augmented body techniques in the modeling of the 
rotational dynamics of the spacecraft and space 
platform. 

1.2 SCOPE 


Study effort was concentrated in five main areas. 

a. Development of a new hybrid ML-LQR approach to 
applying optimal control to the mathematical models 
of the rotational dynamics of prototype flexible 
spacecraft and space platforms. 

b. Improvement in the efficiency of developing three 
axis discrete mass models of the rotational dynamics 
of flexible spacecraft and platforms by utilizing 
augmented body techniques. 

c. Development of three axis discrete mass models of the 
rotational dynamics of a prototype flexible 

pacecraft and of a prototype space platform 
incorporating the flexible spacecraft. 
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d. Application of the hybrid ML-LQR approach to the 
attitude control of models of the rotativial dynamics 
of a prototype flexible spacecraft and the space 
platform. 

e. Digital computer simulation of the rotational 
dynamics models of the prototype flexible spacecraft 
with hybrid ML-LQR attitude control and with LQR 
control for comparison of responses to 
disturbances. 

1.3 GENERAL 

This report is comprised of eight sections. Section 2 
describes the application of a new hybrid ML-r.QR approach 
to effecting attitude control of a single axis three body 
model of the rotational dynamics of a typical flexible 
spacecraft. Section 3 presents the development of a 
three axis five body linearized state variable model of 
the rotational dynamics of che flexible spacecraft 
treated in Section 2., utilizing augmented body 
( barycentric) techniques. Section 4 describes the 
extension of the hybrid ML-LQR approach to the attitude 
control of the three axis five body model and conpares 
the responses of the model with this type of control with 
those obtained from the same model with LQR control 
alone. Section 5 reviews some specific aspects of the 
application of hybrid ML-LQR control to the three axis 
five body model. Section 6 treats the development of a 
three axis ten body model of the rotational dynamics of a 
typical space platform consisting of two spacecraft 
interconnected by a deployable truss described in Ivey 
(1-3). Section 7 presents the extension of the bybnd 
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ML-LQR approach to the attitude control of the three axis 
ten body model of the space platform. Section 9 lists a 
number of conclusions and recommendations. Three 
appendices appear at the end of the report. References 
are cited liberally and are listed at the end of each 
S' ction . 

The original RFQ requested that the International System 
of units (designated as SI) be used in the program and in 
any reporting. Torques, moments, angular momentum, 
moments of inertia and distances, however, are stated in 
English units since this was the system of units used in 
presenting all of the vehicle data in the RFQ and the 
truss data in Ivey (1-3). 
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SECTION 2 


2.0 APPLICATION OF HYBRID MULTILEVEL- LQR ATTITUDE CONTROL TO 

A SINGLE AXIS TORSIONAL MODEL 

In Porcelli (2-1)* an example of the application of 
multiple control in the design of an attitude control 
system for a simple flexible space vehicle comprised of 
three rigid bodies that could rotate about a common axis 
was presented. In Chichester (2-2) an approach to 
applying multilevel techniques to this control problem 
was described. This latter approach vas based upon 
techniques developed in Wismer (2-3) and Chichester 
(2-4), In a later memorandum by Kaczynski (2-5) it was 
shown that single axis attitude control of the discrete 
mass model of Configuration 1 of the Space Construction 
Base could be reduced to the control of the three mass 
single axis torsional model treated in Porcelli' s paper 
and linear quadratic regulator (LQR) techniques were 
applied to efrect this control. 

In the work reported here, an approach to applying 
multilevel control techniques to Porcelli' s model is 
presented that is both more broadly applicable and more 
efficient than the one proposed earlier in Chichester 
(2-2). A combination of Gauss-Seidel multilevel control 
techniques and an extension of LQR techniques is applied 
to this model which contains second order time 
derivatives. The paper upon which this section is based, 


*Numbers enclosed in parentheses refer to references 
listed at the end of this section. 
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Chichester (2-6), ls believed to be the first publication 
of the application of this hybrid approach to the control 
of a flexible vehicle. 


The overall approach followed here in applying the 

combination of multilevel and LQR control m«y be outlined 

as follows: 

a. Express mathematical model of vehicle dynamics in 
(linear) state variable form. 

b. Decompose model into a set of temporarily decoupled 
state equators by defining a suitable set of state 
coordination equations. 

c. Construct correspondingly decomposed performance 
ind ax . 

d. Form decomposed Hamiltonian. 

e. Derive constate equations with associated costate 
coordination equations from necessary conditions for 
minimization of Hamiltonian. 

f. Derive control equations from additional necessary 
conditions for minimization of Hamiltonian. 
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g. Combine state, costate and control equations 
pertaining to the ith riqid body in the model into 
the ith two point boundary value (TPBV) subproblem. 

h. Combine the state coordination and the costate 
coordination equations into an overall coordination 
subproblem. 

i. Assemble the subproblems of the last two steps into a 
two level hierarchy with the coordination subproblem 
in the upper level. 

j. Reduce the solution o£ each TPBV subproblem in the 
lower level o£ the hierarchy to the solution o£ the 
Riccati equation and an auxiliary equation. 

2.1 STATE VARIABLE VEHICLE MODEL 


In Figure 2-1, let: 
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Then the normalized torque equations tor this sinqle axis 
torsional model may be written in state variable form as 
f ol lows : 
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2.2 DECOMPOSED STATE EQUATIONS 

The above state equations may be decoupled temporarily by 
defining the following vector coordination variables. 


^1 

3 » 

-2 

u 
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V 

X 

“1 

Jt 

-2 

31 
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State Coordination 
Equations 


(2-9) 
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With these coord mat ion relationships, the ith state 
equation assumes the tol lowing decomposed form: 


x 

-i 


\i*i 


n u 
1 1-1 


a 

-l 


( t ) ; l - 1 , 


< 2-10 ) 


where: a^ 


A 12— 1 


( 1 1 ) 


a,- A n a" ♦ A^a\ 


< 2-12 ) 


—3* A 32^ 


( 2 - 13 ) 


2 . 3 


DECOMPOSED PERFORMANCE INDEX 


The corresponding decomposed performance index to 
minimize position and rate errors as well as control 
effort may be written m the following form: 


3 

P * X 
l* l 


f 

O 


i, P dt 
t 1 


( 2 - 14 ) 
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Pi 
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- 5di )T 0i ( iU - *dp + 


( 2 - 15 ) 
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i£di * id “ desired angular position and rate lor ail 
three bodies 

2,4 DECOMPOSED HAMILTONIAN 

The specific form ot the Hamiltonian utilized depends 
upon the type oi second level control that is to be 
employed. Feasible, nonleasible and Gauss-Siedel* second 
level control formulations and their respective 
advantages and disadvantages are discussed on pp. 12-16 

*Wismer's name for this approach. 
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of Wismer (2-J). <i*uss-Scidei second level control was 
chosen because it does not require that each subsystem 
have at least as many controls as constraints and it Joes 
not require the use of jradient techniques. From Wiamer, 
pp. 37 - 38 , the decomposed Gauss-iie ide 1 Hamiltonian may be 
written in the following form: 


3 

M - I H < 

l-l l 


T T 

H - (x - x ) Q ( x - x ) + u H u 
1 - i -<11 1 -1 — <3 i -1 1-1 
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T Xu 

+ 3 ( A X + 13 u +• A S + A 3) 

-2 22-2 22-2 21-2 21 -2 


T l T u 
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-2 -2 -1 -2 -2 -3 


14 



U • < X 

3 — 3 
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( 2-21 ) 
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CUSTATK KQUATIONS 

The eostate equations may be written aa tuliows: 
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0. * null vector of appropriate dimension. 
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2.6 


CONTROL EQUATIONS 


The control equations are written in the following form: 


i T 

* 2R U + B \ * 0 (2-27) 

riu i-i ii— i - 

-i 


1 -IT 

u * R B. \ i*l , 2 , 3 (2-28) 

—i 2 l ii— l 


2.7 TPBV FORMULATION OF FIRST LEVEL SUBPROBLEMS 


Substitution of the control equations into the decomposed 
state equations yields the following: 


1 -1 T 
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(2-29) 
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First level subproblem i is formed by associating 
equation (2-23) with equation (2-29). The vector-matrix 
form of this subproblem is expressed as follows: 
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where: 



B . . p . ■ 
ix 1 





Si 


+ 2 Qi*id 


T T 

£i {t o } “ *io 


I0 i-l (t o ) ' °i-l (t o ) ] 


^(tg) * 0 » [0, 0] 


Equation (2-30) and its associated boundary conditions 
constitute a set of three vector two-point boundary value 
( TPBV) subproblems. They may be solved by an extension 
of the linear quadratic regulator ( LQR) techniques that 
appears in Appendix A. 
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GAUSS-SEIDEL SECOND LEVEL COORDINATION SUBPROBLEM 


2.8 


The second level necessary conditions may be written as 
toliows : 

s u * x i*l ,2,3 (2-31) 

— i -l-i 


From 


3 £l 


0 : 


an 


l l 

* 0: s * x l * 1,2,3 (2-32) 

3 v —1 —1+1 


3H 

Is 
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-2 


A T X 
A 12-l 


-3 


A T X 
A 23-2 


(2-33) 
(2-34) 
f 2-35 ) 
(2-36) 


These equations are assembled into the coordination 
subproblem appearing at the apex of the subproblem 
hierarchy depicted m Figure 2-2. 
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FIGURE 2-2: SUB PROBLEM HIERARCHY WITH GAUSS-SEIDEL 

' SECOND LEVEL COORDINATION FOR SINGLE AXIS 
MODEL 
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SECTION 3 


DEVELOPMENT OF A THREE AXIS GAUSS-SEIDEL MULTILEVEL 
ROTATIONAL DYNAMIC MODEL OF A FLEXIBLE SPACECRAFT 

INTRODUCTION 

In the work presented in this section Gauss-Se idel 
multilevel modeling was applied to a three axis five body 
representative of the Space Construction Base (SCB). The 
general approach included augmented body techniques 
presented in Cornell (3-1), (3-2), and Lipski (3-3), 

(3-4) followed by multilevel modeling techniques 
described in the previous section and Wismer (3-5). The 
procedure utilized was as follows. 


1) Approximate the flexible vehicle with rigid bodies 
interconnected by a spring-hinge suspension. 

2) List any symmetry conditions that may apply to the 
elements of the model. 

3) Determine the vector locating the center of mass of 
each rigid body with r«. pect to its barycenter. 

4) Determine the vectors locating each adjacent hinge 
with respect to the barycenter of each rigid body. 

5) Derive augmented moment of inertia matrices 
equivalent to the effective moment of inertia of 
every rigid body in the model with respect to each 
barycenter. 
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6) Write the effective torque balance about each 
barycenter with rigid body angular rates, as 
vector state variables. 

7) Express the equations resulting from 6) in vector- 
matrix form. 

8) Invert the coefficient matrix of w in the model of 7) 
to obtain the state variable rigid body angular rate 
model , 

9) Express time rate of change of relative Euler angles 
in terms of rigid body angular rates to obtain the 
state variable Euler angle model. 

10) Aggregate the state variable Euler angle model with 
the state variable rigid body angular rate model to 
obtain the state variable rotational dynamics 
model . 

11) Decompose tbe state variable rotational dynamics 
model . 

3.1 RIGID BODY-SUSPENSION APPROXIMATION OF FLEXIBLE VEHICLE 

The typical flexible space vehicle depicted in Figure 
3-1, Configuration 1 of the Space Construction Base, was 
approximated by an assembly of five rigid bodies 
interconnected by a spring hinge suspension as shown in 
Figure 3-2. The symmetry relationships for this model 
are as follows. 
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m b “ m 4 * lfl 3 * m 2 ' r 4 ' l i * f 2 C -*—1 > 

£3.3 * - £j.2 £55 - £33 * -Lh “ - J ^2 * ~ £db 

LOCATION OF CENTER OF MASS WITH RESPECT TO BARYCENTEK OF 
EACH RIGID BODY 

Under tnese symmetry conditions tne total mass ot tne 
model is uiven by: 


5 

M » l m * m + 4m ( 3-2 ) 

i«l 1 1 2 


d * 0 d 

-i - -2 


m + m 
I 2 

r = -d 

M -22 -3 


—4 


m I + 3m 2 
M 


r 
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(3-3) 


3.3 


LOCATION OF ADJACENT HINGES WITH RESPECT TO BARYCENTEK OF 
EACH BODY 


—12 “ + —12 ’ —12 


•d 


13 


-24 -2 + -24 


2m^+ bm 2 
_ 


(3-4) 


r„.. ** — d-,. 

-22 -3b 


2b 


3m., 


-22 -2 * Z .22 m HT ^22 " 3 3 


m. 


—44 " ^4 + ^44 * FT £22 ' ~^b 


3,4 GENERATION OF AUGMENTED INERTIA SUBMATRICES 


The etiective moment ot inertia matrix with respect to 
the barycenter ot each rigid body may be expressed as 
toliows ; 


h * 4m 2 D u D u 


h * *2 - m 2 S 2°2 - (m l * 2m 2 )fi 22 l5 22 “ m 2 6 24 S 24 * C 


( 3--j ) 


r 4 _ m 2 D 4 D 4 - <m l f 3m 2 )D 22 D 22 * C 


where: tor _a_ * (a x , ay, 

■v 

a x b * Ab 

f 


a z ) T and J3 * (b x , by, b z ,) T 



a a 0 
Y x 


(3-6) 
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HIGH) HODY ANGULAR RATF EQUATIONS 


i . S 


Equations m terms of tne riqid body rates of angular 
rotation, , may be obtained by equating the rate of 
change of angular momentum of each rigid body comprising 
the model to the sum of torques acting on that body at 
its barycenter. Under the assumption of small angular 
displacements the linearized equat ion for the central 
riqid body in the model may be expressed as follows. 


1 1^1 + Mi5 12 + V * * ^5* 


(3-7) 


-1 f -12 + -13 * ^l-l 


u u (~2 ¥ id ~ i i 




where: 


t. • « vector sum of external torques on ith riqid body, 
tv. a vector sum of torques from iunqe ] on ith rigid body, 

^ J 

f^ * vector sum of external forces on ith riqid body. 


A correspond inq equation was written for each of the 
other riqid bodies in the model. The entire qroup of 
riqid body angular rate equations may be expressed in 
vector state variable form as follows: 
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f l ; j and Aj j * 3 x 3 auhmat rices within symmetric 
matr ic ies , 

The submatrices, A^ , are derived from the original rigid 
body angular rate equations utilizing the symmetry 
properties derived for the hinge location vectors 
relative to the barycanter of each rigid body. 
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3 . b SUSPENSION EQUATIONS 


He suspension equations express the torque acting upon 
each rigid body due to the vector sum ot the torces 
acting on the hinge joints adjacent to the body. They 
may be expressed m vector matrix torm as lollows. 


t H « LG OJ + LK 
— s— a— 


(3-12) 


where: 
h r 
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(3-16) 


3.7 EULER ANGLE RATE EQUATIONS 

The Euler angle rate equations express the first time derivatives of 
the Euler angles of body 1 and the relative Euler angles between 
adjacent bodies of the vehici model in terms of angular rates of 
rotation of the rigid bodies comprising the vehicle. For small angular 
displacements the linearized Euler angle equations are written in the 
following form. 
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X 

a * Klo 


( 3 - 17 ) 


where: K 


~ 0 -I 0 I 0 

0 0 -I 0 I 

-I I 0 0 0 

-I 0 I 0 0 

1 0 0 0 0 


.8 STATE VARIABLE ROTATIONAL DYNAMICS MODEL 


( 3 - 18 ) 


The linearized rigid body angular rate equations and the linearized 
Euler angle rate equations may be combined into a state variable 
rotational dynamic model of the following form. 

A 

ui — GLC u) + GLK a + Gu 
— S — 9— — 

( 3 - 19 ) 

A 

a - Kb 


where _u = _t = (u]r, u^, . . . , = actuator torque vector and the 

remaining vector variables were defined earlier in this paper. If 
these vectors and matrices are expanded in terms of their scalar 
components and elements, respectively, they may be rearranged in such a 
way that scalars associated with a specific axis of the model are 
grouped into the same subvector. Each coefficient matrix appearing in 
the resulting state variable model possesses a high concentration of 
non-zero elements on or near its principal diagonal and few non-zero 
elements distant from this diagonal. This lightly coupled state 
variable model may be written in the following form. 
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x « Z (A, x + B u ) 3 ■ 1/2/3 k « 1,2,3 (3-20) 
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u, a ia tj 1 (k = 1,2,3) have scalar expansions of the 

7S. — is. 

same form as 
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tS2 and ts^ have scalar expansions of the same form as that of _tsj_ with 
the subscripts y and z in the place of x, respectively. 
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is a 5 x 4 matrix ot the same torm as the 5x4 partitioned matrix, 
L, with "1" in the place ot each 

3x3 submatrix, "I". K gi Is a 4 x 5 matrix of the same form as the 4 
x 5 partitioned matrix, Kg, with k]_ 2x the L'^ ace °£ each 3x3 
submatrix, K gl2 , and similar substitutions tor the remaining 
sutmatrices. C gl is a 4 x 5 matrix similarly developed trcm the 4 x 5 
partitioned matrix, C g . The matrices K g2 and C g2 are similarly 
expanded with the subscript, y, in the place ot x as are the matrices 
K g 3 and Cg^ with the subscript, z, in the place ot x. 

Kj^ » 0 for k j. For k * j, it is a 5 x 3 matrix of the same form as 
the 5x5 partitioned matrix, K, with "1" in the place of each 3x3 
subnatrix, "I". 
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MULTILEVEL STATE VARIABLE MODEL 


J.y 


The subscripts appearing m the state variable rotational 
dynamics model, equation (3-20) correspond to the axes ot 
the model. The matrix coetlicients , A^ and 3^, 
represent mteraxial coupling when k j. This model may 
be recast into multilevel form by decomposing it into a 
series of submodels. Decomposition temporarily supresses 
the mteraxial coupling m the overall model producing 
three single axis submodels and a coordination 
submodel. It is effected by writing the following 
coordination equations which constitute the coordination 
subproblem. 


a 

-j 


3 

z 

k-1 


(A d + B s ) 
jk-j jk-j 




■Hk 


j =* 1,2,3 (3-21) 

k*3 » 1,2,3 (3-22) 


Equation set (3-22) is in the form of Gauss-Seidel 
coordination as presented in Wismer (3-5). 

Substitution of equation sets (3-21) and (3-22) into the 
state variable rotational dynamics model of equation 
(3-20) yields three submodels of the following form. 

lj a A jj£j + B j j£j + ( t) j * 1,2,3 (3-23) 
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The four submodels developed in this section may be 
assembled into the two level hierarchy shown in 
Figure 3-3. 

When the relatively light coupling between the axes of 
the above model was neglected, it could be written in the 
following form. 

V + i - 1 ' 2 ' 3 (3 - 24) 
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COORDINATION SUBMODEL 
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FIG. 3-3: SUBMODEL HIERARCHY FOR LINEARIZED DECOMPOSED 
STATE VARIABLE ROTATIONAL MODEL OF VEHICLE 
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SECTION 4 


4.0 COMPARISON BETWEEN HYBRID MULTI LEVEL-LQR AND LQR CONTROL 
APPROACHES APPLIED TO THREE AXIS FIVE BODY MODEL 

4.1 COMPARISON OF SIMULATION RESPONSES 

When the scalar equations comprising the original 
linearized state variable model were rearranged so that 
all equations pertaining to the same axis were grouped 
into the same submodel, the resulting overall model 
contained three submodels with light coupling between 
them. LQR control was applied to a model obtained by 
ignoring the coupling between the submodels associated 
with each of the three axes. Hybrid ML- LQR control was 
applied to a vehicle model that retained the interaxial 
coupling. Both controlled vehicle models were simulated 
on a digital computer. Their responses to an initial 
small angular displacement of one unit about the x axis 
and .1 unit about both the y axis and the z axis were 
plotted , 

Responses of the decoupled model with LQR control appear 
in Figures 4-1 and 4-2. Angular displacement of the 
central body and angular displacement of the outer body 
of the solar panel relative to the body representing the 
inner portion of the solar panel are plotted on common 
axes m both figures. From the plots, it appears that 
the oscillations of the solar panels are reflected in the 
transient responses of the central body. 
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FIG. 4-2: RESPONSES OF THREE AXIS FIVE BODY MODEL WITH LO D CONTROL 
(2 AXIS) 


Response* of the coupled model with hybrid ML-LQH control 
appear m Figures 4-3 through 4-6. In thin series of 
plots the relative displacement of the outer body of the 
solar panel is separated from the plot of the angular 
displacement of the central body to facilitate direct 
comparison between the responses for LQR control and 
hybrid ML-LQR control. From these plots it is evident 
that the two control approaches produce very nearly the 
same responses to the same initial displacements. 

The LQR control approach has been compared with that of 
hybrid ML-LQR control on the basi3 of the responses of 
the digital simulation of the two controlled systems to 
the same set of initial angular displacements. tn the 
sequel, each of these control approaches is summarized 
to facilitate comparison of the two. Each approach is 
applied to one of the models previously presented and 
developed from techniques presented in Lipski (4-1), 
(4-2), Hooker and Margulies (4-3) and Chichester (4-4). 
Comparison of the analyses associated with the two 
control approaches is based in part upon Tiffany (4-6). 
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RESPONSES OF THREE AXIS FIVE BODY MODEL 
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FIG, 4-5; e, RESPONSES OF THREE AXIS FIVE BODY MODEL 
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FIG. 4-6: *1 RESPONSES OF THREE AXIS FIVE BODY MODEL 


4.2 


LQR CONTROL OF THE AXIALLY DECOUPLED THREE AXIS FIVE BODY 
MODEL OF A FLEXIBLE SPACECRAFT. 
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4.2.2 


Generation of Three Hamiltonians and Their Minimization 


Since the axes of the state variable model of the vehicle 
are decoupled from each other, optimal control is applied 
to each axis independently. Optimal linear quadratic 
regulator control of the attitude of the vehicle with 
respect to its jth axis corresponds to minimization of 
the jth Hamiltonian whica incorporates the integrand of 
the quadratic performance index pertaining to control 
about the jth axis. This Hamiltonian may be expressed as 
follows: 

H. - Pj + 4 ( A.x j+ B j u j , (4-2) 


where: Pj 


the integrand of the jth quadratic 
performance index 




_x_j d = specified value of 2Lj 


Qj a* state variable error weighting coefficient 
matrix 

Wj u =* control energy weighting coefficient 
matrix 

_X_j * jth costate vector 
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The ]th costate equation is obtained from a necessary 
condition for minimization of the 3 th Hamiltonian. 






(4-3) 


Another necessary condition for minimization of the 
Hamiltonian yields the control equation. 


4.2.3 



-1 T 
W . B A 

3-3 


(4-4) 


Formulation and Solution of jth TPBV Sub-Problem 


Substitution of the control equation into the state 
equation (4-1) and association of the costate equation 
given by equation (4-3) with the resulting equation 
yields a two point boundary value (TPBV) problem to be 
solved of the following form. 


where : 


x. * A-x.+R.A. 
-J 3-3 3-3 


-1 T 

R. * -B .W ■ B . 

3 3 )U ] 


(4-5) 


x.(t ) x. (initial boundary condition) 
- 30 = -30 

A_.(t f ) = 0 (final boundary condition) 
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Solution of this TPBV problem is effected by assuminq 
that: 


x . 

-0 


K. (t)X. 


(4-6) 


Differentiation of equation (4-6) with respect to time 
and substitution of the resulting expressions for ^ and 
into equation set (4-5) transforms solution of the 
jth TPBV subproblem to solution of the jth Riccati 
equation for K j . 


V - ( Vi + Vi + Wi + V 


1,2,3 


where : 


-1 T 
R ■ * B.W./b! 
1 1 ] 


(4-7) 

(4-8) 


wT^bTk - x • 


F j*i 


K j ( t £ ) * (0) 


(4-9) 


Substitution of equation (4-6) into the control equation 
yields the relationship for obtaining the optimal 
feedback control utilizing the solution to the jth 
Riccati equation. 


4,3 HYBRID ML-LQR CONTROL OF THE THREE AXIS FIVE BODY MODEL 

OF A FLEXIBLE SPACE VEHICLE WITH INTERAXIAL COUPLING. 

4.3.L State Variable Model with Coupling 

The axially coupled state variable rotational dynamic 
model of Space Construction Base Configuration I was 
presented in equation set (3-20). In this set of 
equations, the subscripts identify the axes of the model 
and terms involving both subscripts in which k/j 
represent interaxial coupling. 

4.3.2 Decomposed Model 

Application of Gauss-Seidel multilevel control to this 
system begins with decomposition of the state variable 
model to temporarily suppress this coupling. The 
resulting coordination equations were presented in 
equation sets (3-21) and (3-22). The corresponding 
decomposed model of the rotational dynamics is 
represented by equation set (3-23). Comparison of 
equation sot (3-23) with equation set (3-24) for the 
decoupled model reveals that the decomposed model has 
additional terms on the right hand side due to interaxial 
coupling . 

4.3.3 Generation of Decomposed Hamiltonian and Its Minimization 

Optimal control may be applied to this decomposed state 
variable model by extension of the procedure described 
above for applying LQR control to the model in which the 
three axes are decoupled. 

The modified procedure involves the generation of a 
decomposed performance index and a decomposed Hamiltonian 
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corresponding to the decomposed form of the state 
equations. The application of necessary conditions for 
the minimization of the decomposed Hamiltonian yields a 
set of costate equations and costate coordination 
equations as well as a set of control and control 
coordination equations. 

The state, costate and control coordination equations are 
combined into an overall coordination subproblem. Inputs 
to this subproblem are the state, costate and control 
vectors while its outputs are the various coordination 
vectors . 

4.3.4 Formulation of TPBV Subproblems 

When the control equations are substituted into the 
corresponding state equations and each costate equation 
is associated with its corresponding state equation, the 
resulting set of two-point boundary value subproblems to 
be solved is written in the following form: 


x-» A • • x • + R ■ X • + a.{t) 
-3 3 3-3 3-3 ~3 


— j* Q j*j Aj + bj(t) 


(4-10) 
3 * 1,2,3 


where : 




-1 T 
-B . .W B . . 
33 U3 JJ 


A A 

a.j ( t) , bj(t) * functions of the coordination variables 
Xj(t Q ) » Xj 0 (initial boundary condition) 

Xj(tf) * 0. (final boundary condition) 
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4.3.5 Extended LQR Solution of TPBV Subproblem 

The presence of additional terms on the right hanu side 
of the TPBV problem represented by equation set (4-10; 
compared with the TPBV problem represented by equation 
set (4-b ) requires an extension ot the LQR approach. in 
this extension ot LQR control, it is a-.sumed that tne 
costate variables may be written m the following form. 

If K j(t)x + m^t) (4-11) 

where : 

(t) * time varying matrix tn be determined. 
m^(t) * time varying vector to be determined. 

Differentiating equation (4-11) with respect to time and 
substituting the resulting expressions for X.j an ^ A-j 
equation set (4-lU) transforms the jth TPBV subproblem to 
a pair ot Riccati-type equations to be solved tor Kj and 
m r 

K • “ - (K_ A_ _ + .V^K_ + K_ R_ K_ + Q ) j * 1,2,1 (4-12) 

3 3 33 33 3 D D J j 

i + (K j R : + A jj>£ + K 3 k 2 + bj - 0 (4-13) 

where 

K j (t f ) - [0] 

mj (t f ) « 0 
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Solution of this pair of equations alone} with equation 
(4-11) leads to qeneration of the feedback, control needed 
to minimize the portion of the overall performance index 
associated with the jth axis subject to values of the 
coordination variables obtained from the coordination 
subproblem. Other outputs obtained from the jth TPBV 
subproblem include the state vector/ x^ , and the costate 
vector/ A 

— j . 

4.3.6 Construction of Subproblem Hierarchy 

Due to the relationship between their respective inputs 
and outputs/ the coordination subproblem and the TPBV 
subproblems may be assembled into the subproblem 
hierarchy shown in Figure 4-7. Solution of the overall 
optimal attitude control problem is then attained by 
iteration between the subproblems in the hierarchy. 

4.4 COMPARISON BETWEEN THE TWO APPROACHES 

Had the LQR control techniques been applied directly to 
the coupled model/ each of the coefficient matrices would 
have had dimensions three times as large as they are for 
the coefficient matrices associated with either the 
application of LQR control to the axially decoupled model 
or the application of hybrid ML-LQR techniques to the 
coupled model. Hybrid ML-LQR control/ therefore, reduces 
the problem of optimally controlling the atti'-ude of the 
coupled vehicle model to a form closely approximating the 
control of its attitude with respect to each of its axes 
independently with periodic adjustments by the 
coordination subproblem to account for the effects of 
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mteraxial coupling. Part of the price of attaining thi 
result is greater complexity in the control system and 
the addition of iteration between the subproblems of the 
hierarchy in the generation of the optimal control 
solution. However/ this iterative model of solution 
implies an adaptive capability not present in the 
conventional LQR control approach. 



state, costate and control 


Second Level Coordination Subproblem 
(State, coatate and control coordination equations) 



FIG, 4-7: SUBPROBLEM HIERARCHY FOR HYBRID MULTI LEVEL-LQR 
ATTITUDE CONTROL OF FIVE BODY MODEL 
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SECTION 5 


5.0 REVIEW OF SIMULATION RESULTS OF THREE AXIS FIVE BODY 

MODEL OF A FLEXIBLE SPACECRAFT 

An internal memorandum by N.O. T:'fF a ny (5-1) was reviewed 
and is summarized herein. The purpose of this report was 
to document the results obtained to date With the 
simulations of multilevel and LQR control at the Bendix 
Engineering Development Center in support of this study 
contract and to outline fruitful areas for further 
work. 

Plots of the controlled responses of uncoupled and 
lightly coupled forms of the linearized five body state 
variable rotational model of a flexible space vehicle to 
initial angular displacements were generated on the 
Bendix Engineering Development Center digital computer. 
The prototype vehicle whose attitude was controlled is 
depicted in Figure 3-1 while its corresponding 
topological diagram in which each flexible appendage was 
approximated by two rigid bodies interconnected by a 
spring hinge suspension appears in Figure 3-2. The 
subproblem hierarchy developed for the application of a 
combination of multilevel and an extension of LQR control 
techniques to the coupled form of the model is shown in 
Figure 4-7. The uncoupled linearized five body model was 
derived from the coupled model by ignoring the light 
coupling between its three spatial axes. LQR attitude 
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control was then applied to each axis independently. In 
ti»e application of hybrid ML-LQR attitude control the 
coupled model was partitioned by spatial axes and the 
coupling between these axes was temporarily suppressed as 
explained in Section 3.9. 

Plots of responses of the central body of the vehicle to 
angular displacements about the x, y and z axes appear in 
Figures 5-1, 5-2 and 5-3 respectively. The responses of 
the uncoupled model with LQR attitude control applied 
independently about each axis and of the coupled model 
with hybrid ML-LQR of control to the same initial angular 
displacements are plotted with common sets of 
coordinates. Those responses pertaining to the uncoupled 
model are distinguished by small circles. It is evident 
that the angular responses of the central body of the 
coupled model with hybrid ML-LQR control compare very 
closely with the corresponding responses for the 
uncoupled model with LQR control applied to erch axis 
independently. 

Multilevel control temporarily suppresses interaxial 
coupling effectively reducing the problem of controlling 
the coupled model to the problem of controlling the 
decoupled model with later coordination. 
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Additional response plots were drawn for cases in which 
the control laws generated for the nominal controlled 
system were applied to a system in which the values of 
damping coefficients and the spring constants of the 
model were changed. Plots of the angular attitude 
responses of the central body of the model about its x 
axis appear in Figures 5-4, 5-5 and 5-6 in order of 
increasing perturbation of these parameters. In each 
figure the response of the perturbed system is identified 
by small triangles with the responses of the unperturbed 
system included to facilitate comparison. The principal 
differences between the two classes of responses are 
slightly increased damping and phase shift for the 
responses of the perturbed system. In general, however, 
the control of the vehicle appears quite insensitive to 
such mismatches between the model upon which the 
generation of the control law was based and the system to 
which it was applied. 

5.1 SENSITIVITY OF CONTROL LAW TO SYSTEM PARAMETERS 

The use of a control law depends upon a suitable model of 
the system to be controlled. Suppose that the system to 
be controlled is represented by the state equation 

x * Ax + Bu (5-1) 

while the control law is based on a model of this system 
represented by 

.i it i 

x » A x + B u (5-2) 
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The question arises as to the effect of differences 
between equations (5-i) and (5-2), These may differ in 
many ways including the dimensionality of the system. In 
the case of the five body model, the sensitivity of the 
control to some of the parameters involved was 
investigated. From equation t3-20), the A matrix is of 
the form 


(5-3) 

Here, K s involves spring constants, C g involved damping, 
and G involves mass distribution and geometry. B 
depends upon G only, and the K in equation (5-3) is 
independent of system parameters. 

An investigation was conducted, in which the spring 
constants and damping constants were subjected to 
variations, and the control law as developed without 
those variations was applied. That is, A* in equation 
(5-2) was fixed, and A in equation (5-1) was varied. 

Each damping constant c* contributing to A' was modified 
to c * (1 + DC ) c 1 , and similarly, k « (1 + DK)k' . The 
effects on the central body angle ^ for variations of 
DC, DK * 0,1, 0.5 and 1.0 are shown in Figures 5-4, 5-5 
and 5-6, respectively. The conclusion is that the 
control law is quite insensitive to these parameters, in 
a macro sense. This study was conducted on the 
uncoordinated control only, but the results are 
applicable to the partitioned control also, since the 
coupling is weak. 


A 


GLC GLK 
s s 


K 
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5.2 


CHOICE OF PARTITIONING 


The five body three axis partitioned control converged 
after two iterations, due to the weak coupling between 
axes. The 5 body controL system must operate on 10x10 
matrices, with a resulting penalty in computation time. 
The choice of partitioning was guided by intuition and 
insight. A more systematic approach ls desirable. The 
choice made WLth the five body model attempts to increase 
convergence by satisfying the following independence 

A A 

conditions. If a^ = £ and b ^ 5 _0, then m_j * 0 and the 
jth subsystem is completely controlled by Kj alone, and 
is independent of jc_j . This condition will occur if A ^ j , 
R^j and Q^j are zero for jj*i. Further partitioning may 
be desirable. 

5.3 CENTRAL BODY CONTROL 

In the studies discussed here, control was applied to the 
central body only. Furthermore, for the five body model 
no constraint tending to control flexible appendage 
motion was utilized. The investigation so far has been 
concerned with the practicality and implementation of 
multilevel control via partitioning in general, and has 
not addressed the question of its effectiveness in 
controlling the flexible vehicle. The tools required to 
investigate that aspect of the problem are now developed, 
except for the means for selecting the weights necessary 
to best accomplish the desired control of those states 
representing the flexible behavior. This study assumes 
that all states are available. 
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SECTION 6 


6.0 THREE AXIS TEN BODY MODELING OK THE ROTATIONAL DYNAMU-a 
OF A FLEXIBLE SPACE PLATFORM 

6.1 RIGID BODY-FLEXIBLE SUSPENSION APPROXIMATION OF MODULAR 
SPACE PLATFORM 

The typical platform consisting of two spacecraft 
interconnected by a deployable truss* was approximated by 
an assembly of ten rigid bodies interconnected by a 
spring hinge suspension as shown in the topological 
diagram of Figure 6-1. From this figure/ the following 
symmetry relationships between the rigid bodies of the 
model were written. 


;n 5 - m 4 - m 3 » m 2 £j . 3 - -rj . 2 

r 5 m E 4 * X 3 * 1 2 £55 * £33 * £24 " “£21 * “£35 

Additional symmetry conditions were obtained by comparing 
the truss portion of Figure 6-1 with the description of 
the truss appearing in Ivey (6-1). 

m 8 - m 6 I 8 - I 6 £88 “ £66 * “£89 * “£67 

m 9 M m 7 l 9 m l l £99 * £77 * “£9,10 * “£78 


*A more detailed development of the rigid body spring hinge suspension 
approximation of the deployable truss appears in Appendix B. 
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Selected masses of rigid tiodies compr is mg the model were 
related to each other by ratios as follows* 

m^ « a r ;n 2 ^7 * b r m^ m^g « °r m l 


Then the total mass of the model was express sd In the 
following form 


10 

M • Z m « ( 1 +-2 fc> +c ) in +(4+2a )in 
i*l i r r l r 2 

6.2 LOCATION OF CENTER OF MASS WITH RESPECT TO COK'IECTION 

BARYCENTER OF EACH RIGID BODY 

The connection barycenter of a rigid body may be defined, 
following Hooker and Margulies (6-6), as M the new center 
of mass obtained by 'loading' each joint of the body with 
the residual mass of the system connected to that joint." 
With the aid of the topological diagram of Figure 6-1, 
the symmetry conditions and the ratios between the masses 
defined above, the vectors representing the location of 
the center of mass of each rigid body with respect to its 
barycenter are written as follows. 
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INTERCONNECTED BY DEPLOYABLE TRUSS 




Ijl » -(I L_Uk) r « _ (2b r +c r )m l +2a r m 7 2b „+c 

M 16 w £l6 “ 

1+2 +g 


-16 > 


(l+2b +c )m, + (l+2a_ 


“2 » - W - ^>£ 22 - 


“22 * % > ~£ 2 2 


^ 4-^5 - -a- ^-)r . + t 3 + 2 a r )m ? 


£22 > -£ 2 2 


, r (4b +2c )m + 3a m 

d, » - 1 - r g 1 r 2 

■* !• 5T 1*66 


(l-2b r -c r )m i + (4_ a ^) m 


l-2b-c 

> £ r 

l + 2b r +c r 66 


( 6 - 1 ) 


(3b +2c )m +2a m 
3 - - 1- — r r 1 r 2 

-7 L M 


(1-b -c )m +4m 
r r 1 2 


1-b -c 

s r r 
? — 1 

l+2b +c —77 

r r 


- 1 - 2b r m l + ( 1+ c )nu+(4+a )nu 

M £66 = " “ r 2 


L-*-2b r +c r £66 


73 



2b m 

dg ■ - (I 


( l+o r )m L + (4+2a r )m 2 
- 


-77 


l+c r 

' “ - 11 


Vi 

-10 (1 ~ M ) ~I0,10 


(l+2b r )m^ + (4+2a r )m 2 
- 


-L0,10 


l+2b 


> - 


l+2b r -KJ r -10,10 


6.3 


LOCATION OP ADJACENT HINGES WITH RESPECT TO BARYCENTER OF 
EACH RIGID BODY 


The vectors representing the location of adjacent hinges 
with respect to the connection barycenter of each rigid 
body were expressed in terms of the location vectors with 
respect to the center of mass of the same body as 
follows : 


(2b +c )m,+2a nu 

-12 ^ -1 + -12 = -12 ” -16 


( 2b +c )m, +2a m 
j ^ j . „ _ v . r r 1 r 2 

-13 -1-13 " -13 " M -16 

m 1 +4m 2 1 

-16 = — 1 + — 16 = -16 > l+2b_+c_ -16 


M 


r r 


( 6 - 2 ) 


3m„ 

-22 * - 2 +r 22 = ~ -22 K< 1 
M 33 

3m 

2 

d « d +r = -(2 )r = - d > -2r 

-24 -2 -24 M -22 -35 -22 


m 

2 

d =d+r = — r » - d «1 
-44 -4 -44 M -22 -55 
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d, , * d-+r C£ 

—66 —6 —66 


(4b r +2c r )m L +3a r m 2 2(2b r +c r ) 

r 66 * _ ^66 


M 


l+2b r +c r 


^V7 * V^7 


2m 1+ (8 + a r )m 2 ^ « 2 

^6~^66 * “ M ^66 > I+ib r +c r " ^66 


( 3b +2c )m, +2a nu 3b +2c 

-77 * — 7 + — 77 * ~ -77 v . „ u 71 " -77 


-78 “ — 7 + — 78 3 —7”*— 77 


M ' ' l+2b +c 

r r 

(2+b)m 1 + (8+2a )m,. - / .2+b_) 

E — i — - — r > =■ r__ 

—77 , , nu . - -77 


M 


l+2b +c 
r r 


^sa ‘ 3e + £88 " ^ + £6b 


(2b m.+a m,) 
r l r 2 

M 


2b 


^66 < 


l+2b +c 
r r 


^66 


^89 * ^8 " r -66 


( 2+2b +2c )m. +( 8+3a )m~ -2 ( 1+b +c ) 

^ -** * ± y* S y 

^66 ' ^66 


M 


1+ b+c 
r it 


2b r mi 

^gg * ^9 + ^g * ^e + £.77 M -77 < 

M 


2b. 


l+2b j _+c j _ 


-77 


( 2+2b r +c r )m 1 +( 8+4a r )m 2 

^9,10 3 ^9 + ^-9 , 10 = ^9“-77 3 ” -77 


~2(l+b r +c r ) 

l+2b +c -77 
r r 


C r m l . c r 

-10,10 " -10 + -10,10 + M -10,10 > l+2b r +c r -10,10 


The relative locations of the barycenters, centers of 
mass and hinge locations along the axis of the truss were 
estimated for four sets of values of the mass ratios b r 
and c r . The corresponding sketches, including the 
location vectors, appear in Figures 6-2, 6-3 6-4 and 
6-5. 
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In the first two sketches c r « 1 which corresponds to the 
two spacecraft having approximately the same mass since 
m l >> m 2 ' « 1 corresponds to the total mass of two of 

the bodies attached along the truss being approximately 
equal to the mass of the central body in spacecraft 1. 

The b r * 1/4 condition corresponds to the total mass of 
this pair of bodies being about 25* of that of the 
central body, The same sequence of values of b r was 
repeated in Figures 6-4 and 6-5 for c r » 1/2 which 
corresponds to the mass of spacecraft 2 being about 1/2 
that of spacecraft 1. 

Inspection of Figures 6-2 through 6-5 reveals certain 
common patterns. For both values of c r when b r * 1/4 the 
centers of mass of bodies 6 and 7 remain close to their 
respective connection barycenters. In fact/ for b r * 1/4 
and c r =• 1/2, tne center o£ mass and the barycenter o£ 
body 6 coincide while the center of mass remains very 
close to the connection barycenter for body 7. The 
greater distances between the centers of mass and their 
corresponding barycenters generally occur close to the 
extreme ends of the truss in each case. Decreasing the 
ratio of the mass of the second spacecraft to that of the 
first, c r / reduces the distance between the barycenter 
and the center of mass of the central body of the first 
spacecraft and increases the distance between the 
barycenter and the center of mass of the second 
spacecraft. These effects agree with physical reasoning 
since reduction of this ratio reduces the influence of 
the combined mass of the truss and the second spacecraft 
compared with that of the first spacecraft while it 
increases the influence of the combined mass of the first 
spacecraft and the truss compared with that of the second 
spacecraft. 
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AUGMENTED INERTIA MATRICES 


b . 4 


In the sequel the following form of the vector cross 
product is utilized frequently. Let 


a 


(a x , a y , a z ) T and b 


(b , b , b ) T then a x b * Ab 
x y z 


where s 


A : 






(b-3) 


With the aid ot the symmetry conditions and the sets of 
vectors developed in the preceeding sections and the 
augmented body techniques presented m Lipski (b-4), 

(6-5) the effective moment ot inertia matrix with respect 
to the connection barycenter ot each rigid body m the 
model is expressed in terms of the above listed mass 
ratios as follows: 


r * 

n 


Z 1 " Wl " 4m 2 D l2 D 12 “ 


(2b c+ c ]f )m 1+ 2a r ^ 


V “ r 2 - ™2°2°2 -( < 1+2b r + = r >V <2+2 V m 2 )fi 22 E 22 "“AAl ’ X 3* 

I 4 * - I 4 - m 2 D 4 D 4 -{ (l+2b r+ o r )m 1 +(3+2a r ) m2 )d 44 d 44 - l b * 

(6-4) 

V * X b -WA -("> 1 +4» 2 )5 66 3 66 - ( (2b r +0 1; )m i +a r m ;! 1 
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I 7 * - I 7 - ~ * (4M r Jm^ ,D 77 D 77 

- 1 - <‘> r «-V )m l* a r m 2 - B 78 5 78 

‘a* * ! a * V2V9 * ' <l*V m l H 4 M r , " l 2 D «Au 

- ,b r w r ,m l B 89 5 B9 


V * [ 9 ■ b r m l D 9 D 9 ( lrt V m l f (4+ia r )m 2 I u 99 D 99->-' r m l D 9 L0 D 9,10 
l 10* ’ I I0 “ C r m L B 10°10 


Under the assumption that the mass of body 1 is much 
larger than the mass of each rig id body in the model of 
the solar panels, >> m 2 and the differences between 
the rotational moments of inertia for each rigid body are 
bounded as follows: 
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S3 



In earlier sections the relative locations of the 
baryeenters, centers ot mass and hinge Locations along 
the axis of the truss were estimated for tour sets of 
values of the mass ratios b r and e r . In the present 
3ecr ion the differences between the rotational inertia 
matrices with respect to the center of mass and with 
respect to the connection barycenter for each rigid body 
in the model are obtained for the same four sets of 
values of b r and c r , The normalised differences between 
the inertia matrices for the rigid bodies along the axis 
of the truss appear in Table 6-1. Since R^ j may differ 
considerably from for j*i, comparisons between rows 
along a given column have more validity thAn do 
comparisons between columns of a given row of the 
table , 


84 



TABLE b-L 

NORMALIZED DIFFERENCES BETWEEN ROTATIONAL * NEKT I A MATRICES 

FOR BODIES ALONG TRUSS AXIS 


I - I * 
b X b 

I 7 - 1/ 

1 V 

I,- l>/ 

I - I * 

1 10 no 

Vbb*6b 

m L R 77 R 77 

m L R bb R bb 

m l l5 7/ S T7 

m I R L0,I0 R lQ / 10 








6.5 


RIGID BODY ANGULAR PATE EQUATIONS 


Equations in terms of the rigid body angular rates of 
rotation, wu , were obtained by equating the rate of 
change of angular momentum of each rigid body in the 
model to the vector sum of torques acting on that body at 
its connection barycenter. Under the assumption of small 
angular displacements, the symmetry conditions, the 
vectors developed in equation sets (6-1) and (6-2) and 
the moment of inertia matrices developed in equation set 
(6-4), the rigid body angular rate equations were written 
as follows: 


« . ■ • • • XU 

Iw+MlDDw+DDw+DDoj+DDto+D £ ( D u) ) J 

1-1 12 22—2 13 33-3 12 44-4 13 55-5 16 j-6 jj-j 


HHH~ ~ ~ ~ ~ 10 ~ 

» t + t + t + t + D f - D f - D f-D f-D Z f-D f 

-1 -12 -13 -16 1-1 13-3 12-4 13-5 16 j-6 -j 12-2 


* , ~ 10 
Io)+M(DDuj + D Z (D to ) + D D to J 

2-2 22 12-1 22 j»3 jj-j 24 44-4 

3*4 


H H ~ 10 

■t+t+t+Df-D Z f - D f 
-2 -22 -24 2-2 22 j*l -j 24-4 

j 2,4 


* . ~ 10 ~ ~ . 

Iu + M( D D a) + D Z (D co ) + D D to 1 

3-3 33 13-1 33 jj-j 35 55-5 

3-2 

j-3,5 
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H H ~ ~ 

t + t + t + D £ -D 
-3 -33 —35 3-3 33 


1U 

Z 


£ -D f 


1 -j 35-5 


PS, 


( 6 - 6 ) 


I u + M( D D jo + D D u + D Z (D u ) J 

4—4 44 12— L 44 24-2 44 j j— j ^ 

3*3 

3 3*4 

H ~ - .10 

* t + t + D £ + D Z t 

—4 —44 4—4 44 — j 

3*1 


3*4 


X u + Mf D D w + D D u + D 
5-5 55 13-1 55 35-3 55 


10 ~ . 

,Z (Du) 
3*2 33-3 

3*3,5 


1 

j 


H ~ ~ 10 

■t + t + D £ - D If 

-5 -55 5-5 55 -3 

3*1 

3*5 




. ~ 5 ~ i-1 

Xu + MfDDu + D Z (D -a ) + D Z {D u ) 

1— i 11 16-1 ii j j — j ii 3/3+1— j 

j*2 j*6 


10 

+D Z (Du)) 

1/1+ 1 . ^ j j— j 

3*1+1 


H H - i-1 

*t+t + t +D £ - D Z (£ ) 

— i — i,i — i,i+l i— i ii j«l — j 


D 

i,i+l 


10 

s 

j*i+l 


(t.) 

-3 


1 * 6, 7, 8, 9 
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* . , ~ 5 ~ . - 9 - , 

I w + M D D -j + D Z (D <d ) + D Z (D at ) 
10-1.0 10,10 16-1 10,10 jj-j 10,10 j, 5+1-3 

3-2 j -6 


H ~ ~ 9 

t+t +D1-D E f 

-10 - 10,10 10-10 10,10 -j 

j-1 


10 

where M * Z m * m +4m +2m +2m + m 
i-l i 1 2 6 7 10 


(l+2b +c )m +(4+2a )m 
r r ' r 2 


6.6 VECTOR STATE VARIABLE FORM OF RIGID BODY ANGULAR RATE 

EQUATIONS 

The rigid body angular rate equations of equation set 
(6-6) were recast into the following vector matrix state 
variable form. 


Aw * t + t H + Bt + 3i * A _1 (t + t H > + A" i Bf 


where ; 


. , . T . T . T J 

oj ** ( co a ^ ^ f • • • / ^10 


(6-7) 
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^ ■ angular rotation rate ot ith rigid body 

T T T T 

~ * ( —X ' —2 r * * * ' —X 0 ^ 

£ » net torce acting on ith rigid oody 


t 


ii 


H 



t 


T 

2 ' ’ * * ' 



) 


T 


vector sum ot control and disturbance 
torques on ltd rigid body. 


,.H « T , . H > T , , h H . T , H .T 

12 X3 + 16 } ' 22 + — 24 


H .T , H T , H .T , H T . H ,T +( H % T 

(£33) + 3 5 ' ( ^44 } ' —55 ' '—66 -—67 ' 


35' 


‘-44 


—55 ' —66 


,,H . T ,.H , T .T., h ,T . . H .T ,T 

(t 77 ) +(t 78 ) , {ty Q ) +(tg 9 ) , ( Jt 9 9 ) + ( Jig ,10 r —1 0 r 1 U 


H 

_t“ * torque ot 3 th hinge on ith rigid body. 


A 11 A l, 10 

• • 

• • 

A * . . * 10 x 10 partitioned 

. . coetticient matrix 


A 10,l • • * * A 10,10 
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B 11 B 1 , 10 

• • 

• • 

B « . . ■ 10 x 11) partitioned 

. . coetticient matrix 


b lo,i • * • ,B 1U,1UJ 

B^ are 3x3 submatrices 


b.7 EXPANSION OF COEFFICIENT SUBMATRICES IN TERMS OF 

AUGMENTED MOMENT OF INERTIA MATRICES AND LOCATION VECTORS 
WITH RESPECT TO CONNECTION BARYCENTERS 


The following matrices within the coefficient matrices, A 
and B, were expanded with the aid of the Lett hand sides 
ot the rigid body angular rate equations listed in 
equation set (6-6). 


A. • * I • l * i, 2 1 • • • 1 10 

11 X t 


/V* 


A 12 - MD i2 D 22 “ A 21 

, T 


rsj v 


A 14 * MD I2 D 44 * A 41 


*vj rv 


A U - MD 13 D 3J - A 31 


A 15 ' MD U D bb - A 51 


A lj =“ MS 16 S jj - A jl' 3-6-7,8,9,10 


A 23 

* M °22 S 33 

T 

* A 32 

A 24 * 

MS 24 B 44 

,T 

A 42 

A 2j 

* M5 22 S jj 

T 

a* A . 

J2' 

j»5,6,7,8,9,10 



A 34 

“ M °33 S 44 

T 

- A « 

A 35 

MD 35 D 55 

« a t 
A 53 
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A 3j « MD 3 3 D ^ ^ * A? 3 , 3-6,7,8,9,10 

A 4j " m 6 44 6 3 j * A j 4 ' 3-5,6,7,8,9,10 (b 

A 5j “ M °55 6 3j " A 35' 3-6,7,8,9,10 

A L j » 1-6,7 ,8,9, j-i+l, l+2, ... ,10 


B ^ —D ^ , 1— 1,2,... ,10 


B 12" _D i2‘ B U B 13* ‘ D 13 * B 19 B lj" -°j.„i3-6,7,8,9,10 


B 2L“ ” D 22* B 2J B 24* -D 24 B 2j“ " D 22 ! 3-5,6,7,8,9,10 


B 31* _D 33 = B 32“ B 34 B 35* -D 35 B 3j’ “ D 33 ; 7-6,7,8,9,10 


B 4l“ " D 44 ' B 42“ B 43 B 4j * “ D 44 J 3-5,6,7,8,9,10 


B 51 “ -B 55“ B 52* B 53 = B S4 B 5j “ ~ B 55 ! 3-6,7,8,9,10 

B ij" - S ii ! i-6 ,7,8,9,10 Bij » -D lji+1 ; 1-6, 7, 8, 9 

j— l,2,.».,i— 1 j— l+l,i+2,...,10 
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SUSPENSION EQUATIONS 


b.8 


The suspension equations represent the torque acting upon 
each rigid body due to the vector sum of the forces 
applied at the hinge joints adjacent to the body. They 
were expressed in vector matrix form as follows. 


t H * LC j + LKi 

— S— 3 — 


where: 


H.T H.T H.T T 

( ' 9 • •• 9 (t l(J ) 


(b-LO) 


(b-11) 


££ » vector sun of hinge torques acting on ith rigid body. 


T T T T 

r t * * * f —]_0 ^ 


£ * rate of rotation ot ith rigid body. 

T 

3 ( 'jJ a 'I) * ijj ) 

-l v lx 7 iy' iz' 


T T T T T T T T T T ,T 
- “ ( -l 7 -12 7 -13 7 -24 7 -35 7 —16 7 -^67 7 -^78 7 ^89 7 -^9,10 ) 


( V V V 


^l 7 ^l 7 ^l * mertially referenced Euler attitude angles of body 1. 

T 

"X .* ( A<t / AB / Aili ) 

-rj v F ij' 13 7 v rj' 


/ A0 ij 7 * relative Euler angles ot body j with 
respect to body l. 

Coefficient matrices, L, C s and K s are expressed in terms 
of 3 x 3 submatrices in Appendix C. 
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6.10 


EULER ANGLE RATE EQUATIONS 

The Euler angle rate equations express the tirst time 
derivatives ot the mertially referenced Euler attitude 
anqles ot body 1 and the relative Euler attitude angles 
between adjacent bodies of the model in terms of angular 
rates of rotation of the bodies comprising the model. 

For small angular displacements, the linearized Euler 
angle equations were written m the toLiowmg torm.. 


Kuj 


( 6 - 12 ) 


where: 


■I * 



T T T T T T T T .T 
-13 ' -24 ' -35 ' -16 ' 1^7 ' 1*78 ' ^89' ^9,10 J 





STATE VARIABLE ROTATIONAL DYNAMICS MODEL 

The linearized rigid body angular rate equations, 
linearized suspension equations and linearized Euler 
angle rate equations were combined into a state variable 
rotational dynamics model of the following form. 
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( »' - 1 3 ) 


L * GLC *> + GLK +■ Gu 

— s~ a— — 

• 

x • Ku 
where : 

* rigid body angular rate vector, 

a * Euler angle vector, 
u « actuator torque vector, 

and the matrices G, b, K, C s and K R are expanded in terms 
of 3 x 1 submatrices in Appendix C. 

When these vectors and matrices were expanded in terms of 
their scalar components and elements, respectively, they 
were rearranged in such a way that scalars associated 
with a specific axis of the model were grouped into the 
same subvector. Each coefficient matrix appearing in the 
resulting state variable model was written in the 
following form. 

3 

Xj - ]C E 1 (A jk x k + B jk u k ) j- 1 , 2,3 k » 1 , 2,3 ( 6 - 14 ) 

where: 

m m m 

x k - («£/ a k ) ( 6 - 15 ) 
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3 * 


1 *]k'’k sk j ^jkA^ak 


<] ]k 

v _ T _r ° r \ 

]k 

101 

m - 


101 - 10 x 10 
z».*ro irwatnx 


-1 * '"lx' "dx' 


lOx 


)‘ 


-2 " U ly' “V' UJ 10y ) 


T 


-3 * U lz' U 2Z , " ^, "IOZ 5 


T 


u. (k » 1,2,3) has scalar expansion ot the same form as -u, , 

— K 


ii 


A* 


12' 


A* 


l 3 ' 


A* 


A* 


r l 


35 


A$ 


lb' 


A$ 


b7 f ^7U' 






~2 


( V Ao 12' 




13' jy 24' 35 ' 


4o 


lb' 


At, 67' ia 78' 


89' 


• iw 9 / io ) 


( + 1 / -* 4 # 


1' * v 12' >H 13' ^24' 




35' 


-V, 


lb' A V7' "* v 78 ' “*89 ' “*9,10 






I’suiq the definitions ot 0/ L, K, C g , and K g as given in 
Appendix 0% is a 10 x 9 matrix of the same form as the 
10 x 9 partitioned matrix, L, witn "l" in the place of 
each submatrix "I" and "0" in place of each subnatrix 
" 10]". K g j_ is a 9 x 10 matrix of the same form as the 9 
x 10 partitioned matrix, K g , with "k^x in the P lace o£ 
submatrix "K^" and similar substitutions for the 
remaining submatrices. C gJ _ is a 9 x 10 matrix similarly 
developed from the 9 x 10 partitioned matrix C g , The 
matrices Kg^ are similarly expanded with the subscript, 
y, in the place of x as are the submatrices K g3 and Cg 3 
with the subscript, z, in place of x. * [ u 1 for k /* 

j. For k * j, it is a 10 x 10 matrix of the same form as 
the partitioned matrix K, with "1" m the place of each 
"I" and "0" in the place of each "[0J M . 
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SECTION 7 

7.0 HYBRID Ml LTlLEVEL-LtfR ATTITUDE CONTROL OF THE ROTATIONAL 
DYNAMICS MODEL OK THE SPACE PLATFORM 

7.1 DECOMPOSED STATE VARIABLE MODEL 

The subscripts appearing in the state variable rotational 
dynamics model, equation (6-14), correspond to the axes 
ot the model. The matrix coefficients, A^. and B, k , 
represent interaxial coupling with k * j . This mod**! was 
recast into multilevel form by decomposing it into a 
series of submodels. Decomposition temporarily 
suppresses the interaxial coupling m the overall model 
producing three smqle axis submodels and a coordination 
submodel. It is effected by writinq the following 
coordination equations which constitute the coordination 
subproblem . 

3 k k 

a * 2 (Ad + 13 3 ) ] * 1 , ., 3 (7-1) 

-3 k-i jk-] ]k-] 

kM 

k. k 

d * x 3 * u k * j * 1,2,3 ( 7-2) 

~3 -] -k 

Equation set (7-2) is in the form of Gauss-Beide 1 
coordination as presented in Wismer (7-L). 

Substitu’ on of equation sets (7-1) and (7-2) into the 
state variable rotational dynamics model of equation 
(6-14) yields three submodels of the following form: 
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J • i,2,\ 


The four submodels developed in this 'section were 
assembled into tne two level hierarchy shown in 
Figure 7-1. 

Corresponding to the decomposed state variaole equations 
of this: model, a decomposed performance index and a 
decomposed Hamiltonian were constructed. 

7.2 DECOMPOSED PERFORMANCE INDEX AND HAM I ETONIAN 


The decomposed performance index for the application of 
optimal attitude control to this vehicle may be written 
as follows: 

t (7-4) 

J f 

j - r. , r dt 

t i 


where s 


1 > 

J 


T T 

^ X , -X ) 0 (x - X ) i* V2 U W u . 

-J -] d 2 -] -J d -J ]U-] 



pcespecif ied desired value of x^ 


q. » symmetric positive semidefinite state variable error 
- 1 weighting coefficient matrix. 

W. * symmetric positive definite control energy weighting 
^ coefficient matrix. 
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COORDINATION SUBMODEL 
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FIG. 7-1: SUBMODEL HIERARCHY FOR LINEARIZED DECOMPOSED 
STATE VARIABLE ROTATIONAL MODEL OF VEHICLE 





The corresponding decomposed Hamiltonian expressing the 
energy state of the system is written as follow 


H « E H. 
j -1 3 


(7-5) 


where : 


H « P + X 

j j 


T r 

A x 

•j L jj-i 


3 

k k 1 

+ B u +• Z { A d + B s) 
jj-j . . jk-j jk-j J 

k*l 


k* j 


(7-6) 


3 r k T j kT j I 

S, (p.) (x,-d J ) + (v.) (u.-s ) 
*1 L -3 -3 -k -3 -3 k J 


k» 

k* j 


_\j * jth vector costate variable (Lagrange multiplier) 





and 


Gauss-Seidel type vector coordination 
variables 


7.3 


COSTATE EQUATIONS 


The costate equations result from application of the 
following necessary condition for minimization of the 
decomposed Hamiltonian: 


X * 

-j 


3H 

3x 


T 

A X - Q (x - x ) + b (t) 
DD-j 3 “j -jd -3 


(7-7) 


where : 
b. (t) 


3 

E 

k=l 

k*j 


k 
p . 
-3 


(7-8) 

X (t ) = 0 (final 
— j f — 

boundary conditions) 
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7.4 


CONTROL EQUATIONS 


7.5 


The control equations result from application of the 
following necessary condition for minimization of the 
Hamiltonian? 


3H 

3u 

-j 


0 






3 

+ E 

k*l 

k*j 


k 

v.) 

-3 


(7-9) 


COMBINATION OF jTH STATE, CONTROL, AND COSTATE EQUATIONS 
INTO THE jTH TPBV SUBPROBLEM 


Substitution of the jth control equation into the jth 
state equation yields a new state equation in terms of _Xj 
and Xj. Association of the jth costate equation with the 
equation representing the combination of the state 
equation and the control equation results in the jth two- 
point boundary value (TPBV) subproblem. 


x. * A. .x. + R. X ■ + a. (t) 
-3 33-3 3-3 -3 


X.* - Q.x. - A?.X • + b. (t) 
-3 3-3 33-3 -3 


where: 


(7-10) 


3j(t) 


3 

* E 

k*l 

k*j 


(A., d^ + B., a^ - B. . \fc) b.(t) 
3k-3 jk— 3 33-3 -3 


°* d “k»l 
k*j 


k 


R. 


-B. .WTW. 
33 3U 33 


x.(t ) * x. (initial boundary conditions) 
—3 o —jo 

_X. (t f ) =0 (final boundary conditions) 

3 t 
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Solution of the subprobleme requires an extension of the 
LQR techniques to reflect the presence of the additional 
terms on the right hand side of the corresponding TPBV 
subproblems. This extension, which was presented in 
Chichester (7-2), beginu with generalizing the assumed 
form of in terms of _Xj to the following. 

i,j(t) * K^(t)Xj(t) + m j v •’ ) 

wh^re K j ( t ) and m_j(t) are a time variable matrix and 
vector, respectively, to be determined. The assumed Corm 
of results in the generation of two simultaneous 
Riccatl type equations to be solved for Kj(t) and 
nuj(t). Substi tut: on of these solutions into equation 
(7-11) yields the feedback control to minimize the 
performance index at time tf. It should be noted that 
the corresponding solutions to the TPBV subproblems are 
obtained subject to the value of the coordination 
variables d^, s ^ , p_ ^ , and v^. 

7.6 ASSEMBLY OF COORDINATION EQUATIONS INTO THE COORDINATION 

SUB PROBLEM 


The remaining necessary conditions for minimization of 
the Hamiltonian result in additional coordination 
equations to be satisfied. 



d* 



0 





(7-12) 


3H n ^ k _T , 

^ ° \j>k 


3H 

3_Sj 


- v 3 ? 
— D 


B kj— k 
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The top two of the above coordination equations are 
associated with the state equations while the bottom two 
are associated with the costate equations. All of the 
coordination equations are assembled into an overall 
coordination subproblem. Inputs for this coordination 
subproblem are state vectors, x_j , control vectors, _u_.j , 
and costate vectors, _Aj c , obtained as solutions of the jth 
and kth TPBV of the hierarchy respectively. 

7.7 CONSTRUCTION OF THE SUBPROBLEM HIERARCHY 

The original problem of generating the optimal attitude 
control contours for the prototype flexible space vehicle 
has been transformed into a series of interrelated 
subproblems to be solved. Due to these interconnections 
it is evident that these subproblems can be assembled 
into a hierarchy with the coordination subproblem at its 
apex. This hierarchy is depicted in Figure 7-2. 

Solution of the overall optimal attitude control problem 
is attained by iteration between the levels of the 
hierarchy. The modification of the space platform's 
configuration by adding or removing modules is treated at 
the first level of the hierarchy as presented in 
Chichester 
(7-2) . 
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Second Level Coordination Subproblem 
(State, costate and control coordination equations) 
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First Level TPBV 
Subproblem 1 
(State, costate 
and control eq 'ns) 


First Level TPBV 
Subproblem 3 
(State, costate 
and control eq'ns) 


FIG, 7-2: SUBPROBLEM HIERARCHV FOR HYBRID MULTI LEV EL- LQR 

ATTITUDE CONTROL OF TEN BODY MODEL 
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SECTION 8 


8.0 CONCLUSIONS AND RECOMMENDATIONS 

During FY'81 a new hybrid approach to the application of 
multilevel control techniques was developed and utilized 
in the attitude control of a series of prototypical 
models of flexible space vehicles. This hybrid approach 
consists of a combination of the techniques of multilevel 
hierarchical (ML) control and an extension of linear 
quadratic regulator (LQR) techniques to treat the 
modified two point boundary value subproblems resulting 
from the application of multilevel techniques. 

This combination of ML and extended LQR techniques was 
first applied to a single axis model of the rotational 
dynamics of a prototype flexi: le spacecraft. The 
resulting model of the spacecraft and its control system 
was simulated on a digital computer. 

The same combination of techniques was then extended to 
the attitude control of a three axis five body model of 
the rotational dynamics of the same prototype flexible 
spacecraft. This extension entailed the approximation of 
the flexible spacecraft by five rigid bodies 
interconnected by a spring hinge suspension and the 
development of the model of the rotational dynamics in 
terms of connection barycenters. The resulting model was 
linearized and then recast into a state variable form 
amenable to the application of time domain control 
methods such as multilevel and LQR approaches or a 
combination of the two. When each of the vector-matrix 


106 



equations comprising this model was expanded into its 
scalar components a regrouping of these scaler eouationn 
according to spatial axis was found to yield a model in 
which interaxial coupling was light. 

Two approaches were followed in applying attitude control 
to the above model. Under the first approach the light 
coupling between the axes was eliminated and LQR control 
was applied to each axis. Under the second approach the 
interaxial coupling was temporarily suppressed by 
applying multilevel control techniques in conjunction 
with an extended form of LQR techniques. The three axis 
five body model of the rotational dynamics of the 
flexible spacecraft and its attitude control system wa3 
simulated for each of the two approaches described to 
facilitate comparison. 

The hybrid ML-LQR approach was further extended to the 
attitude control of a three axis ten body model of the 
rotational dynamics of a prototype space platform 
comprised of two spacecraft interconnected by a 
deployable truss. This model was derived utilizing the 
same techniques as those used in the development of the 
three axis five body model described above. 

As was the case with the five body model the ten body 
model was linearized, cast into state variable form and 
expanded in terms of scalar components. The scalar 
equations of this model were then rearranged in such a 
way that all equations associated with a given spatial 
axis were placed in the same group. The hybrid ML-LQR 
approach was then applied to the attitude control of this 
model. Simulation of the model with hybrid ML-LQR 
attitude control is pending at this time. 
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8.1 


CONCLUSIONS 


The following conclusions are based upon single axis and 
three axis modeling, control and simulation ot a 
prototype flexible spacecraft and modeling and control 
formulation for a three axis representation of a 
prototype space platf irm. 

a. Multilevel techniques reduce a mathematical model of 
high dimension (involving large numbers of variables 
and equations) to a hierarchy of submodels each of 
which is of smaller dimension except for the overall 
coordination submodel at the top of the hierarchy 
which generally is of much simpler form than the 
original model. 

b. Multilevel techniques effect a similar reduction of 
a control problem of high dimension to a hierarchy 
of subproblems to be solved. 

c. Since application of multilevel techniques to a 
control problem containing second order differential 
equations becomes Ineffective if an attempt is made 
to reduce any of the first level subproblems to a 
dimension less than two, such application leads to a 
vector-matrix formulation of the subprcblems in the 
bottom level of the hierarchy. 

d. Application of multilevel techniques to an optimal 
control problem decomposes the overall two point 
boundary value (TPBV) problem to be solved to a two 
level hierarchy of subproblems with TPBV subproblems 
on the lower level and a coordination subproblem on 
the upper level. 
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e, Each of the TPBV subprobleins on the lower level of 
the hierarchy resultin'] from the application of 
multilevel techniques to an optimal control problem 
contains coordination variables on its right hand 
side requirtnq an extension of linear quadratic 
regulator ( LQR) techniques for its solution. 

f, Tlu* combination of multilevel and extended LOR 
techniques described above represents a new hybrid 
approach to the application of optimal attitude 
control to a mathematical model of rotational 
dynamics of a spacecraft. 

q, Diqttal computer simulation of a single axis model 
of the rotational dynamics of a prototype flexible 
spacecraft with a hybrid ML-LQR attitude uontroL 
system produced very nearly the same responses as a 
simulation of the same single axis model with an LQR 
attitude control system. 

h. Development of a linearized three axis five body 
state variable model of the rota ional dynamics of 
the prototype flexible spacecraft was rendered more 
efficient by utilization of augmented body 
techniques. 

i. Expansion of the three axis five body modeL in terms 
of scalar variables- >nd equations revealed that 
grouping these equations according to spatial axes 
yielded a model with a high degree of block diagonal 
dominance in its coefficient matrices with weak 
inheraxial coupling. 
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j. Digital computer simulation ot the three axis five 
body model with hybrid ML-LOR attitude control 
produced responses closely similar to those 
resulting from computer simulation of the same model 
with axes decoupled and L.OR control applied to each 
ax is, 

k, Development ot a three axis ten body linearized 
state variable model of the rotational dynamics ot a 
prototype space platform consisting of two 
spacecraft interconnected by a deployable truss lead 
to approximation of the truss by four rigid bodies 
serially connected by a spring hinge suspension. 

l, Application of hybrid ML-L.QK techniques to the 
attitude control of the throe axis ten body model of 
the space platform yielded a subproblem hierarchy of 
a form similar to that developed tor the three axis 
five body model of the prototype flexible 
spacecraft . 

m. The hybrid ML-LQR approach to attitude control of 
flexible spacecraft is an original concept that 
appears quite promising for application to a oroad 
variety of spacecraft rotational dynamics models. 
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Kh’COMMfcNDATIUNS 


The t'ol lowing items* are recommended tor tuture study 

concerning control of flexible spacecraft. 

a. The tnree axis ten body model of the rotational 
dynamics of the prototype space platform with hybrid 
■nul 1 1 level-LOK attitude control should i« simulated 
on a digital jomputer m order to evaluate its 
responses to a variety of standard perturbations. 

To this end the NASTRAN data for the MSKC/hybrid 
truss supplied by NASA should be utilized to 
generate a dynamically equivalent rigid body spring 
hinge suspension model of the truss comprised of 
tour serially connected rigid bodies, 

b. The relative cont~'l authority available with 
actuators on both spacecraft in the space platform 
should be comparer! with that which is available when 
control actuators are restricted to one of the two 
spacecraf t , 

c. Different schemes for decomposing the original 
attitude control problem into a hierarchy of 
subprobleins should be compared with respect to 
computational efficiency. 

i. Techniques for unprovinq the computational 
efficiency of solving the modified two point 
boundary value subproblem resulting from the 
application of the hybrid ML-LQR approach to the 
attitude control of rotational models of a variety 
of spacecraft should be developed. 
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e. Effects of modeling inaccuracies upon the 
performance of the controlled system should be 
evaluated . 

f. More systematic methods for determining the 
weivjhtlmj matrices utilized in the decomposed 
performance matrices associated with the hybrid ML- 
LQR approach to attitude control should be sought. 
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APPENDIX A 


EXTENSION OF LINEAR QUADRATIC REGULATOR TECHNIQUES TO 
FIRST LEVEL TPBV SUBPROBLEMS* 

From Equation (2-30) 


ii * A u*r 7 ' 1 1 ! 

i - 1,2,3 

h - -2Q j x i -A^ii + b i (t) 


(ai; 


with given x^(t o ) and j^(t f ) 


(A2) 


Assume: 

-i * K i^ t) -i + JSi ( fc ) 


X_- * K . x • +K . x • +m . 
—l l—i l—i —l 


(A3) 

(A4) 


Substituting equation (A3) in equation (Al) 


(A u- k (A5 > 


Substituting equation (A4) in equation (A2) 




T 


-(A^K^O^^-m. (A6) 


*This is similar to an approach presented in Sage (2-7) 
pp. 96-97. 
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Substituting equation (A5) in (A6): 

<V K 1 A U - I K i«i' lK l + A iT K i +2Q l>*i 
1 * -1 A . 

+ ^i^~i~ "J R i — + — i i * £ iA7) 


For solutions for K^(t) and mj_(t) independent of 


Ki«.Aii - £ K.R.- 1 K i+ Aj i K. + 2Q.. [0] 


— i + iK.Ri^m.-K.ai-^ - 0 


{ A8 ) 
(A9) 


Since jL^tj) “ _0_r for solutions independent of _Xj_, 
K^(t f ) » [0] and m^(t^) * £ (A10) 


Also, 

u i ( t > 



~\l -i (t) 




T 


B i i ( K i— i 


+ Sl ) 


(All) 
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APPENDIX B 


SUCCESSIVE APPROXIMATIONS OF MSFC/HYBRID DEPLOYABLE TRUSS 

A perspective drawing of the undeformed MSFC/Hybrid 
deployable truss that appeared in Ivey (6-1) is 
reproduced in Figure B-l. The following specifications 
were received for the truss. The truss is comprised of 
ten segments each of tfhich is a 55-inch cube. Each of 
these segments weighs 16.8 pounds. If the segments are 
numbered consecutively from the left hand end of the 
drawing, four 8000 pound rigid bodies are rigidly 
connected in such a way that they are centered above and 
below the fifth and tenth cubic segments with their 
centers of mass 89 inches above and below the truss as 
shown in Figure B-l. Additional data resulting from a 
vibrational analysis of the truss were presented in Ivey 
( 6 - 1 ) . 
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First Approximation of MSFC Truss 


L, Each cubic segment of the truss was approximated by a 
cube of the same dimensions with the total mass cf the 
segment uniformly distributed within the cube. 

2. The cubic segments were associated into two pairs of 
identical modules as shown in Figure B-2. 

3. One of the module pairs has external masses associated 
with it as depicted in Figures B-l, B-2 and B-4 . 

4. Equivalent principal moments of inertia were obtained for 
each module pair. (Figures B-3 and B-4). 

5. Modules were serially connected by a spring-hinge 
suspension with spring and damping coefficients 
associated with each three degree-of-f reedom hinge shown 
in Figure B-2. 
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FIG, B-3: TRUSS MODULE 1 

TOTAL MASS = = 4(168) CENTERED AT PT, 0 

IOg 


4= 27,5 IN, 

THE PRINCIPAL MOMENTS OF INERTIA ARE AS FOLLOWS : 



I Y = 68 miLj 

T 

I z = 68 

y 
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Principal Moments of Inertia 


I x = ,rn 2 (1 l + 1 2 )Z + ^l 1 ! = m 2 [ ( 1 1 * 1 2^ * 
l z I 111 ! 1 ! 




B-120 


Second Approximation of MSFC Truss 


1. Truss may be approximated by four serially connected 
rigid bodies as depicted in Figure B-5. 

2. Each rigid body of the second approximation has the same 
principal moments of inertia as the corresponding truss 
module in the first approximation of the truss depicted 
in Figure B-2. 

3. The second approximation has the same three degree-of 
freedom hinges/ spring constants and damping coefficients 
as the first approximation. 
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APPENDIX C 


DEFINITIONS FOR MATRICES APPEARING IN THE TEN BODY STATE VARIABLE 

ROTATIONAL DYNAMICS MODEL 


L 


91,30 


1.930, L * * * * 930, 30J 


,-1 


G U G L2 G U 

G 21 G 22 G 3 3 

G 31 G 32 G 33 


- Cllk' 
G 2k - <i2k' 
G 3 k ■ ( a 3k ' 
2.ji * ( 9ji' 

A 2 " <«j2' 

2 j 3 ’ (g j 3 ' 


24k' 

2 7 k 

2sk' 

2 8 k 

26k' 

29 k 

9 j4 ' 

9 j 7 

9 j 5 ' 

9 j 8 

9 j 6 ' 

9 j 9 


* * * ' 223 k 5 

* **' 2 . 29 k 5 

* * * ' 230 k 5 

* * * ' 9 j28 5 
• • • / 9 j 29 ^ 
***' 9 j 3 Q 5 
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DAMPING COEFFICIENT MATRIX 
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